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water  was  dripped  down  the  surface  of  a  0. 5A in. -thick,  /LNg-chilled  SOFI 
panel.  A  preliminary  study  of  solar  melting  produced  little  rundown  or 
ice.  Frost  densities  of  less  than  0.1  g/cnw  (6  lb/ft^)  were  obtained, 
except  in  the  test  in  which  water  was  dripped  onto  the  surface. 

Validation  by  the  present  experiments  of  predictions  based  on  a  heat 
balance  model  applies  only  to  natural  convection  in  the  laminar  flow 
regime.  A  large  fraction  of  the  flow  on  the  space  shuttle  fuel  tank  will 
be  in  the  turbulent  regime  because  of  the  height  of  the  tank.  Larger- 
scale  tests  would  be  required  to  validate  the  turbulent  heat  transfer 
coefficients. 
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I.  INTRODUCTION 


Operation  of  the  Space  Transportation  System  (STS)  from  Vandenberg  Air 
Force  Base  (VAFB),  rather  than  from  the  warmer  Kennedy  Space  Center  (KSC),  has 
raised  the  question  of  whether  ice  will  form  on  the  foam-insulated  external 
tank  containing  cryogenic  propellants.  Ice  accumulation  and,  in  particular, 
falling  ice  loosened  by  the  launch  represent  a  potential  hazard  to  the 
orbiter;  a  launch  will  be  postponed  for  certain  predetermined  ice  condi¬ 
tions.  The  existing  mathematical  model^  of  heat  transfer  and  ice  formation 
predicts  that  launches  would  be  delayed  more  often  at  VAFB  than  at  KSC. 

The  purpose  of  this  study  was  to  obtain  subscale  experimental  data  to 
validate  some  portions  of  the  mathematical  model.  In  addition,  the  present 
study  served  as  a  "pathfinder"  for  larger-scale  experiments. 

Predictions  by  the  mathematical  model  indicate  that,  for  certain  air 
temperature  and  humidity  conditions  that  depend  on  the  wind  velocity,  cloud 
cover,  time  of  day,  and  insulation  thickness,  the  surface  temperature  of  the 
external  propellant  tank  can  fall  below  the  freezing  point  of  water.  Freezing 
occurs  most  readily  when  the  air  is  still  (no  forced  convection)  and  only 
natural  convection  and  radiation  warm  the  surface  of  the  tank.  For  thiB 
reason,  we  decided  to  test  a  simulated  section  of  the  external  propellant  tank 
(ET)  in  a  controlled  environmental  chamber  under  the  still  conditions  of 
natural  convection.  We  controlled  temperature  and  humidity,  and  measured  the 
buildup  of  frost  mass  and  thickness  at  intervals  during  tests  lasting  up  to 
6  hr.  Surface  temperatures  were  also  monitored.  The  measured  frost  accretion 
(accumulation)  rates  and  surface  temperatures  were  compared  with  the  predic¬ 
tions  of  a  simplified  heat  balance  equation. 


II.  APPARATUS 


The  tests  were  conducted  In  an  environmental  chamber  measuring 
6  x  6  x  11  ft  that  was  lined  on  the  inside  with  foam  Insulation.  The  chamber 
has  a  built-in  air  conditioner/blower  in  the  ceiling.  However,  when  the  air 
conditioner  was  operated  for  several  hours  under  high-humidity  conditions,  the 
coils  would  frost  over,  reducing  the  ability  of  the  air  conditioner  to  cool 
the  chamber.  Consequently,  we  placed  two  large-area  heat  exchangers  in  the 
chamber  to  permit  cooling  without  excessive  moisture  extraction  from  the  air. 
The  heat  exchangers  were  cooled  with  liquid  nitrogen  (LN2)-chilled  water- 
glycol  mixtures  to  temperatures  a  few  degrees  above  freezing,  i.e. ,  near  32  to 
35°F.  The  humidity  in  the  chamber  could  be  controlled  independently  with 
three  residential  humidifiers  to  nearly  100  percent  relative  humidity  (RH). 

Test  panels  were  constructed  to  simulate  a  section  of  the  external 
tank.  They  were  made  by  attaching  a  36  x  16-in.  test  sample  of  sprayed-on 
foam  insulation  (SOFI)  to  one  side  of  a  1-in. -thick  aluminum  frame.  The  SOFI 
test  samples  were  obtained  from  Martin-Michoud  and  were  CPR  488  sprayed  onto 
the  1/8-in.  aluminum  sheet.  The  back  side  of  the  frame  was  covered  with  a 
sheet  of  1/8  x  36  x  16-in.  aluminum.  The  aluminum  sheets  were  fastened  to  the 
frame  with  machine  screws  and  sealed  with  RTV  silicone  cement.  The  test 
panels  were  designed  to  hold  LN2,  and  Styrofoam  insulation  was  placed  on  all 
exposed  aluminum  surfaces.  The  SOFI  samples  had  an  original  thickness  of 
between  1-1/2  and  1-3/4  in.  We  shaved  one  sample  to  0.5-in.  thickness  for  the 
first  test  series  (tests  7-17)  and  another  sample  to  0.9  in.  for  the  second 
series  (tests  18-24).  The  SOFI  was  shaved  to  obtain  faster  frosting  rates  and 
to  facilitate  obtaining  frost  samples  from  the  smoother  surface. 

The  panel  was  filled  with  LN2  through  a  l/2-in.  insulated  line  from  a 
commercial  LN2  tank.  Nitrogen  gas  boll  off  was  exhausted  through  a  line  to 
the  outside  of  the  building.  The  pressure  in  the  panel  was  monitored  with  a 
Heise  Bourdon  gage,  and  the  LN2  flow  rate  was  controlled  to  prevent  a  pressure 
greater  than  about  1  pslg.  A  glass  dipstick  containing  two  thermocouples  was 
inserted  through  the  top  of  the  panel  to  detect  LN2  at  1  ft  off  the  bottom  and 
3  in.  from  the  top. 
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Thermocouples  (type  T,  Cu-Constantan)  were  attached  with  Scotch  tape  to 
the  outside  surface  of  the  SOFI,  6  In.  from  the  top  and  12  in.  from  the 
bottom.  The  temperature  profile  through  the  SOFI  was  measured  with  thermo¬ 
couples  embedded  at  various  depths.  In  order  to  place  the  thermocouples  at 
various  depths ,  we  cut  a  triangular  trench  through  the  SOFI ,  its  apex  at  the 
aluminum  plate.  The  thermocouples  were  placed  along  the  side  wall  of  the 
trench  at  three  levels  in  the  0.5-in.  SOFI  and  at  four  levels  in  the  0.9-in. 
SOFI.  The  thermocouple  lead  wires  were  routed  for  several  inches  along  an 
approximately  isothermal  path  to  avoid  temperature  perturbations  by  heat  con¬ 
duction  along  the  wires.  After  the  lead  wires  were  positioned,  the  triangular 
piece  of  foam  was  replaced  in  the  trench  and  sealed  along  the  edges  with  a 
thin  application  of  RTV  cement. 

The  test  panel  was  mounted  vertically  in  the  environmental  chamber,  its 
bottom  edge  approximately  6  in.  from  the  floor.  The  ambient  air  was  sampled 
for  dew-point  temperature  (T^)  measurements  at  four  locations,  two  in  front 
of  the  panel  and  two  in  back.  Those  locations  were  at  the  same  height  as  the 
center  of  the  panel  and  12  in.  from  the  panel  surface.  Folyflo  tubing  carried 
the  sampled  air  to  a  General  Eastern  dew-point  instrument.  Thermocouples  were 
attached  near  the  ends  of  the  sampling  tubes  to  measure  the  ambient  tempera- 

ture  ^ Tamb  ^  * 
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III.  PROCEDURE 


The  Initial  cooling  down  and  filling  of  the  panel  required  approximately 
40  min,  which  is  comparable  to  the  fill  time  expected  at  the  launch  site. 

Most  of  the  thermocouples  were  monitored  with  a  computerized  monitoring  and 
recording  system.  A  few  of  the  temperatures,  including  the  dew-point  tem¬ 
peratures,  were  read  and  recorded  manually.  During  the  cool-down  period,  the 
flow  of  LN2  was  carefully  controlled  to  avoid  any  large  pressure  buildup  in 
the  panel  that  could  have  broken  the  seals.  Once  the  panel  was  nearly  full, 
the  thermocouples  on  the  dipstick  approached  the  temperature  of  LN2,  and  the 
pressure  in  the  panel  oscillated  at  about  5  cps.  While  the  panel  was  being 
filled,  the  Internal  temperature  of  the  environmental  chamber  was  chilled  by 
the  built-in  air  conditioner,  as  well  as  by  the  large-area  heat  exchangers, 
toward  the  desired  operating  conditions.  The  air  conditioner  was  shut  off 
once  the  panel  was  full. 

During  the  course  of  a  test,  the  desired  temperatures  were  maintained  by 
making  small  adjustments  to  the  LN2  flow  rates  through  the  panel  and  through 
the  water-glycol  bath  of  the  heat  exchanger.  Periodically,  the  panel  was  in¬ 
spected  to  determine  the  extent  of  frost  formation.  The  frost  thickness  was 
measured  with  a  metal  scale  scribed  in  0.01-ln.  divisions.  Samples  of  frost 
were  scraped,  with  a  razor  blade,  from  a  3-to-5-in.  x  1.5-in.  area  and  were 
weighed  on  a  triple-beam  balance  with  a  precision  of  10  mg.  The  dimensions 
of  the  sampled  area  were  measured  and  recorded. 

The  above  scraping  technique  worked  well  for  the  shaved  SOFI  panels  but 
might  prove  difficult  on  the  rough  texture  of  an  as-sprayed  panel.  Therefore, 
we  tested  a  technique  in  which  the  soft  frost  was  swept  from  an  area  into  a 
collection  dish  with  an  artist-type  brush.  The  frost  quantities  measured  in 
that  way  were  slightly  less  (10-20X)  than  those  obtained  with  the  razor-blade 
scraping.  The  brush  technique  may  prove  useful  for  as-sprayed  surfaces  on 
which  scraping  would  be  difficult. 
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IV.  RESULTS 


A.  PRELIMINARY  TESTS  1  THROUGH  6 

Preliminary  tests  1  through  6  were  performed  with  panels  sprayed  in  our 
laboratory  with  a  foam  obtained  from  the  Aerospace  shipping  department.  Those 
tests  demonstrated  that  the  large-area  heat  exchangers  were  required  to  main¬ 
tain  the  desired  constant-temperature  and  -humidity  conditions. 

B.  SHAVED  PANEL— 0.5  in.  THICK 

Tests  7  through  17  were  conducted  with  a  NASA-Martin  Marietta-supplied 

SOFI-covered  panel  that  was  shaved  to  a  thickness  of  0.5  ±  0.03  in.  The  set 

of  data  for  test  12,  presented  in  Fig.  1,  indicates  that  both  the  weight  and 

thickness  of  the  frost  increased  linearly  with  time,  starting  45  min  after  the 

panel  was  completely  full.  The  test  lasted  6-1/2  hr  from  the  start  of  the  LN£ 

fill  to  the  last  measurement.  The  slope  of  the  data  gives  a  frost  mass  accre- 

—3  2 

tion  rate  of  9.0  x  10  g/cm  -hr  and  a  thickness  growth  rate  of  0.16  cm/hr. 

The  density  was  fairly  constant  during  the  test  at  0.055  g/cm  as  evidenced  by 
the  linear  plots  of  the  frost  mass  accretion  and  the  thickness  versus  time  in 
the  figure.  The  ratio  of  the  frost  mass  accretion  to  the  thickness  determines 
the  density. 

Figure  2  shows  the  temperature  of  the  SOFI  surface  as  determined  by  a 
thermocouple  held  in  place  by  a  thin  strip  of  Scotch  tape.  (The  tape  has  a 
high  conductivity  compared  with  the  SOFI  and  should  not  perturb  the  measure¬ 
ment  greatly.)  Once  the  LN2  fill  was  started,  the  temperatures  dropped 
rapidly  and  reached  -7  to  -8°C  (19-18°F)  about  the  time  the  Ll^  level  indi¬ 
cator  showed  the  panel  to  be  full.  The  temperature  held  at  that  level  for 
about  45  min  (comparable  to  the  time  delay  before  the  frost  appearance  in 
Fig.  1).  After  the  45-min  delay  time,  frost  became  visible  and  the  SOFI 
surface  temperature  began  to  decrease,  eventually  reaching  -22°C  (-8°F)  by  the 
end  of  the  test.  The  surface  temperature  of  the  frost  was  not  measured; 
therefore,  the  temperature  drop  across  the  frost  cannot  be  determined. 

However,  the  SOFI  temperature  undoubtedly  decreased  because  of  the  insulating 
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effect  of  the  frost.  The  frost  deposition  data  obtained  in  tests  8  through  13 
are  listed  in  Table  1. 

Several  special  tests  were  conducted.  Test  7  was  a  shakedown  test; 
test  14  was  a  water-rundown  test;  and  tests  15  and  16  were  demonstrations.  In 
test  14,  conducted  at  ambient  conditions  of  15°C  and  30  percent  RH,  water  at 
ambient  temperature  was  dripped  onto  the  front  of  the  L^-filled  panel  at  four 
equally  spaced  positions  about  4  in.  down  from  the  top.  Hypodermic  needles 
were  used  to  restrict  the  flow  to  about  1  drop/sec  or  0.7  g/min  (1.5  *  10-^ 
lb/rain)  at  each  position.  The  water  ran  down  the  panel  in  small  streams  about 
3  mm  (1/8  in.)  wide  and  began  freezing  on  the  lower  part  of  the  panel  within 
30  min  after  the  fill  was  completed.  The  ice  thickness  grew  almost  linearly 
with  time  over  the  course  of  the  5-hr  test  to  a  final  thickness  of  about 
1.0  cm  (0.4  in.),  and  a  width  of  approximately  2.0  cm.  The  ice  had  an  approx¬ 
imately  semicircular  cross  section.  The  surface  temperature  of  the  ice-free 
portion  of  the  panel  varied  between  -7°C  (19°F)  toward  the  top  and  -11°C 
(12°F)  toward  the  bottom.  Toward  the  end  of  the  test,  the  surface  area  of  the 
ice  was  sufficiently  large  to  freeze  all  of  the  running  water.  Photographs  of 
the  panel  with  the  four  icicles  are  presented  in  Fig.  3. 

C.  SHAVED  PANEL— 0.9  in.  THICK 

Tests  18  through  24  were  performed  with  the  NASA-Martin  Marietta-supplied 
S0FI  panel  machined  to  a  thickness  of  0.9  in.  Frost  formed  more  slowly  on 
this  panel  than  on  the  0. 5-in. -thick  panel.  The  frost  grew  linearly  in  mass 
and  depth  after  about  the  same  delay  time  (45  min)  as  observed  for  frost 
formation  on  the  0.5-in.  panel;  however,  the  measured  densities  tended  to  be 
higher  toward  the  beginning  of  a  test  than  toward  the  end.  Changes  in  density 
as  large  as  a  factor  of  2  were  observed,  and  most  occurred  during  the  first 
1-1/2  hr. 

Ambient  conditions,  frost  growth  rates,  densities,  and  surface  tem¬ 
peratures  for  tests  18  through  22  are  listed  in  Table  2.  Ambient  temperatures 
ranged  from  4.6°C  (40°F)  to  19°C  (66°F).  Most  of  the  tests  were  conducted 
between  40  and  50°F  and  at  97  percent  RH.  All  of  the  frost  densities  were 
<0.1  g/cnr*  or  <6  lb/ft^. 
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TABLE  1.  EXPERIMENTAL  DATA  FOR  0.5-in.  PANEL 


Test 

Tamb 

Tdp 

RH 

Frost 

Growth  Rate 

Density 

Tsurf 

(°P> 

(°C) 

(°F) 

(*) 

<g/c«2 

-hr) 

(c«/hr) 

(g/ca3) 

(°C) 

(  *F) 

8 

64 

18 

53 

66 

6.8 

X 

10-3 

0.088 

0.080 

-8 

18 

9 

61 

16 

60 

100 

9.2 

X 

l 

o 

— * 

0.150 

0.060 

-3 

28 

10 

48 

9 

48 

100 

4.9 

X 

10"3 

0. 160 

0.031 

-12 

10 

Ua 

59 

15 

57 

95 

14 

X 

10"3 

0.140 

0.090 

-1 

30 

12b 

57 

14 

57 

100 

9.0 

X 

r*i 

1 

O 

0.160 

0.055 

-7 

1 * 

13 

55 

13 

34 

41 

2.7 

X 

ID'3 

0.080 

0.034 

-13 

9 

aConducted  with  fan  blowing  at  aedlua  speed. 
^Conducted  with  fan  blowing  at  low  speed. 
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D.  EFFECT  OF  SIMULATED  SOLAR  RADIATION 

The  effect  of  simulated  solar  radiation  was  examined  in  tests  23  and  24, 
using  the  0.9-in.  panel.  As  depicted  in  Fig.  4,  a  100-W  light  bulb  mounted  in 
a  ref lector/ shield  was  placed  about  7  in.  from  the  middle  of  the  panel.  Some 
of  the  power  to  the  bulb  was  lost  as  heat,  and  only  part  of  the  100  W  was 
radiated.  However,  if  100  W  are  assumed  to  radiate  in  all  directions,  then 

the  flux  on  the  panel  directly  in  front  of  the  bulb  was  ~0.025  W/cm  or  about 

2 

one-fourth  the  maximum  solar  flux  of  0.1  W/cm  . 

Test  23  was  run  at  a  Tam^  of  9°C  (48®F)  and  a  T^  of  45  ±  1.5°F.  The 

*  o 

lamp  was  turned  on  after  a  thin  layer  (~3  ram)  of  frost  (~2.1  *  10  J  g/cm  )  had 

formed.  The  frost  melted  in  about  2  min  and  appeared  to  soak  into  the  shaved 

—3  2 

S0FI  surface.  The  accumulated  frost  was  less  than  the  14  x  10  g/cm 

2 

required  for  water  runoff. 

Test  24  was  a  repeat  of  test  23  except  that  more  time  was  allowed  for  the 

_2 

frost  to  accumulate  before  the  lamp  was  turned  on.  The  frost  (~1.3  x  10 

o 

g/cm  )  directly  in  front  of  the  lamp  melted  within  7  min,  and  that  in  an  area 
about  4  in.  in  diameter  melted  during  the  next  10  min.  Some  of  the  frost 
above  the  lamp  melted  because  of  convection  currents  induced  by  the  hot  lamp 
housing.  A  few  drops  of  water  moved  slowly  down  the  panel,  but  most  of  the 
water  appeared  to  stay  in  place  or  diffuse  into  the  unmelted  frost. 


18 


V.  HEAT  BALANCE  MODEL  PREDICTIONS 


The  three  important  variables  of  surface  temperature  (Tsurf)»  [(frost 
mass  accumulation)/area] ,  and  thickness  were  measured  in  these  simulated-ET 
experiments.  A  heat  balance  model  relates  the  first  two  variables  to  the 
environmental  conditions  and  physical  characteristics  of  the  SOFI.  The 
equations  are  given  in  the  Appendix,  but  a  few  remarks  about  the  input 
parameters  are  appropriate  here.  The  heat  conduction  through  the  SOFI  depends 
on  SOFI  thermal  conductivity;  we  have  used  the  NASA-Martin  Marietta  values. ^ 
The  equation  for  laminar,  convective  heat  transfer  to  a  vertical  flat  plate 

O 

was  taken  from  McAdams.  The  mass  diffusion  coefficient  for  H20  in  air  was 
derived  from  the  heat  transfer  coefficient  by  the  Schmidt  analogy.  The 
radiative  transfer  term  requires  the  emissivlties  of  the  SOFI,  the  frost,  and 
the  foam  insulation  lining  the  inside  walls  of  the  environmental  chamber.  The 
emissivlty  of  frost^  is  close  to  unity  in  the  important  infrared  portion  of 
the  spectrum  that  dominates  radiative  transfer  near  room  temperatures.  The 
foam  lining  and  the  shaved  SOFI  can  be  expected  to  have  large  emissivlties 
because  of  their  small  cellular  structure.  Even  the  unshaved  SOFI  with  no 
open  cells  has  a  high  emissivity,  reported  to  be  0.92.^  We  used  a  value  of 
1.0  in  the  heat  balance  calculations;  using  a  smaller  value  would  have  yielded 
a  slightly  lower  SOFI/frost  Tsurf  and  a  slightly  larger  rate  of  frost 
accumulation. 

The  results  of  the  calculations  for  tests  8  through  13  (0.5-in.  SOFI 
sample)  are  compared  with  the  measured  values  in  Table  3.  Tests  11  and  12 
were  conducted  with  the  air  conditioning  fan  at  medium  and  low  speeds, 
respectively  (no  air  conditioning).  The  measured  frost  accumulation  was 
increased  somewhat  by  that  added  convection.  The  calculated  Tsur£  are 
compared  with  those  measured  in  the  lower  half  of  the  SOFI  after  the  panel  was 
filled  with  LN2  but  before  the  frost  began  to  form.  The  average  ratios  of  the 
calculated  values  to  the  measured  values  for  tests  8,  9,  10,  and  13  were  1.13 
and  1.02  for  the  frost  growth  rate  and  (Tanb  -  Tsurf)»  respectively.  This 
agreement  is  good,  considering  the  uncertainties  in  the  experimental  data  and 
input  parameters,  as  well  as  the  theoretical  assumptions. 
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TABLE  3.  COMPARISON  OF  EXPERIMENTAL  AND  CALCULATED  VALUES  OF  FROST 

GROWTH  RATES  AND  SURFACE  TEMPERATURES  FOR  TESTS  8  THROUGH  13 
(0.5-in.  SOFI  PANEL) 


Test 

Frost 

Growth  Rate 

(mg/cm  -hr) 

(T 

—  T 

amb  surf 

)  CC) 

( meas ) 

(calc) 

(calc/ meas) 

(meas) 

(calc) 

(calc/meas) 

8 

6.8 

6.84 

1.01 

26 

22.6 

0.87 

9 

9.2 

9.7 

1.05 

19 

21.3 

1.12 

10 

4.9 

6.7 

1.37 

21 

23.0 

1.10 

lla 

14 

— 

— 

16 

— 

— 

12a 

9.0 

8.3 

0.92 

21 

21.7 

1.03 

13 

2.7 

2.96 

1.10 

26 

26.1 

1.00 

Avg. 

Runs  8 ,  9 . 

10,  and  13 

1.13 

1.02 

^ests  affected  by  blower. 


The  measured  and  predicted  data  for  tests  18  through  22  are  presented  for 
comparison  in  Table  4.  The  predicted  values  of  (Taml)  -  Tgurf)  fall  very  close 
to  the  measured  values  and  are  well  within  the  uncertainty  of  the  measured 
values.  The  predicted  values  of  the  frost  growth  rates  averaged  about  12  to 
13  percent  higher  than  the  measured  values,  although  the  scatter  of  the 
experimental  data  is  somewhat  large. 

For  the  typical  case,  the  radiation  term  in  the  heat  balance  equation  is 
somewhat  larger  than  the  contributions  of  convective  heat  transfer  and  heat  of 
condensation/fusion.  However,  all  three  terms  are  significant  in  the  calcula¬ 
tion  of  Tsurf«  The  frost  condensation/fusion  rate  depends  directly  on  the 
convective  heat  transfer  coefficient  assumed  by  the  Schmidt  analogy.  Both  are 
only  weakly  dependent  on  the  difference  (Taa^,  -  Tgurf). 
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TABLE  4.  COMPARISON  OF  EXPERIMENTAL  AND  CALCULATED  VALUES  OF  FROST 

GROWTH  RATES  AND  SURFACE  TEMPERATURES  FOR  TESTS  18  THROUGH  22 
(0.9-in.  SOFI  PANEL) 


Test 

n 

Frost  Growth  Rate  (mg /cm  -hr) 

(T 

aab  ”  Tsurf^  ^  c-* 

(meas) 

(calc) 

(calc/aeas) 

(meas) 

(calc) 

(calc/meas) 

18 

— 

— 

— 

16 

15.4 

0.96 

19A 

4  *  1 

4.1 

1.0 

14 

13.7 

0.98 

19 

4.6 

4.4 

0.96 

14 

13.3 

0.95 

20 

3.2 

3.5 

1.10 

13 

14.0 

1.08 

21 

2.6 

3.6 

1.39 

15 

13.9 

0.93 

22 

2.9 

3.4 

1.17 

15 

14.1 

0.94 

Avg  1.12 


Avg  0.97 


VI.  DISCUSSION 


The  test  results  obtained  with  the  0.5-  and  0.9-in. -thick  SOFI  panels 
indicate  that  fairly  accurate  surface  temperatures  and  frost  accumulation 
rates  can  be  predicted  with  a  heat  balance  model  calculation.  The  model  does 
not  predict  frost  density  or  thickness,  however.  Measured  and  calculated 
frost  accumulation  rates  are  plotted  in  Fig.  5  and  agree  well  except  for  the 
two  points  taken  with  the  fan  in  operation.  Although  the  number  of  measure¬ 
ments  is  small,  the  present  data  suggest  that  the  calculations  may  overesti¬ 
mate  the  frost  accumulation  rate  by  about  12  percent.  One  source  of  error  in 
the  measurements  may  be  that  the  convective  flow  down  the  panel  was  influenced 
by  the  floor,  which  was  only  6  in.  below  the  bottom  of  the  panel.  The  cooled 
air  had  to  turn  and  flow  across  the  floor  to  make  room  for  more  air  to  flow 
across  the  SOFI  surface.  Uncertainties  are  introduced  into  the  calculations 
by  inaccuracies  in  the  convective  heat  transfer  coefficient  recommended  by 
McAdams,^  the  mass  transport  coefficient  obtained  by  the  Schmidt  analogy  from 
the  heat  transfer  coefficient,  the  thermal  conductivity  values  for  the  SOFI, 
and  the  SOFI  thickness. 

The  most  important  variable  affecting  frost  accumulation  in  the  various 
tests  is  the  humidity  or  the  partial  pressure  of  f^O  in  the  air.  The  measured 
frost  accumulation  rates  are  plotted  in  Fig.  6  versus  (Famb  -  P3urf)»  the 
difference  in  the  partial  pressure  of  water  in  the  ambient  air  and  the  vapor 
pressure  of  water  at  the  surface  of  the  SOFI  (or  frost).  The  mass  transfer 
coefficient  depends  weakly  on  AT,  the  temperature  difference  between  the 
ambient  air  and  the  surface,  but  only  to  the  one-fourth  power  of  that 
quantity. 

Investigators^  have  predicted  with  some  success  the  density  of  frost 
formation  on  uninsulated  flat  plates.  Huffmann^  at  the  University  of  Dayton 
Research  Institute  (UDRI)  has  described  the  frost  structure  with  a  mass 
diffusion  model  and  applied  that  theoretical  model  to  the  present  problem  of 
frost  formation  on  an  Insulated  plate.  Unfortunately,  the  initial  calcula¬ 
tions  yielded  frost  densities  between  0.2  and  0.5  g/cra^,  whereas  all  densities 
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CALCULATED  FROST  ACCUMULATION  RATE 


1 


MEASURED  FROST  ACCUMULATION  RATE 

•  /  2  . 
m  ,  g/cm  -hr 

Fig.  5.  Calculated  frost  accumulation  rates  versus  measured  frost 
accumulation  rates,  o  “  0.5-in.  SOFI;  A  *  0.9-in.  SOFI. 
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MEASURED  FROST  ACCUMULATION  RATE,  g/cm  -hr 


Fig 


6.  Measured  frost  accunulation  rates  versus  water-vapor  partial 
pressure  (P^  -  P„urf). 
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are  plotted  versus  the  calculated  Tgurf  In  Fig.  7.  They  can  be  fitted  with  a 
straight  line  in  the  range  of  -1  to  -148C  with  a  maximum  deviation  of  ±15 
percent  for  any  point.  They  correlated  somewhat  better  with  the  calculated 
TSurf  than  with  the  measured  ones.  Extrapolations  outside  the  present 
temperature  range  would  have  a  high  degree  of  uncertainty. 

In  this  study,  the  frost  mass  and  thickness  grew  linearly  with  time.  In 
studies  of  uninsulated  vertical  plates,  Nakamura  (see  Dietenberger  et  al.®), 

O  Q 

Cremers  and  Mehra,  and  Kennedy  and  Goodman’  have  observed  frost  thickness 
growing  as  the  square  root  of  time,  a  phenomenon  also  observed  in  forced 
convection  over  uninsulated  tubes.1®  The  growth  rate  decreases  with  time  as 
the  insulating  layer  of  frost  accumulates  on  the  otherwise  uninsulated  sur¬ 
face.  The  insulation  serves  to  raise  the  frost  surface  temperature  and  slow 
the  rate  of  water  condensation/fusion.  In  the  present  case,  the  heat  transfer 
rate  through  the  SOFI  remains  essentially  constant  since  the  insulating  effect 
of  the  frost  is  small  compared  with  that  of  the  SOFI. 

The  delay  time  observed  for  the  start  of  frost  formation  was  about  45  min 
for  both  thicknesses  of  SOFI.  The  various  stages  of  crystal  growth  have  been 
described  in  detail  by  Hayashi  et  al.;11  also  see  Cremers  and  Mehra.®  The 
first  sign  of  frost  was  the  appearance  of  small,  needlelike  crystals  growing 
out  of  the  surface.  They  began  to  branch  out  and  link  together,  gradually 
forming  a  more  uniform  structure.  The  frost  surface  has  a  much  larger  effec¬ 
tive  surface  area,  which  may  facilitate  frost  growth.  The  initial  limiting 
process  to  frost  formation  may  be  crystal  growth  on  the  smooth  surface,  with 
diffusion  of  water  through  the  boundary  layer  being  the  limiting  process  only 
at  longer  times. 

The  running-water  tests  resulted  in  large  icicles  attached  very  firmly  to 
the  SOFI  surface.  They  could  not  be  dislodged  by  hand,  and  no  attempt  was 
made  to  pry  them  off  the  fragile  SOFI.  The  first  drops  of  water  on  the  panel 
moved  somewhat  erratically  down  the  surface  because  they  did  not  wet  the  sur¬ 
face.  However,  once  a  path  was  formed,  succeeding  drops  followed  it.  Once 
the  column  of  water  had  begun  to  freeze,  the  drops  of  water  ran  down  the  ice, 
leaving  a  thin  film  of  water  to  add  to  the  icicle  growth.  An  estimate  of  the 
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ice  thickness  growth  rate  with  the  heat  balance  model  agreed  fairly  well  with 
the  measurements,  representing  the  heat-transfer-limited  case,  in  which 
sufficient  water  is  assumed  to  be  present  at  all  times. 

The  results  of  the  simulated  solar  experiments  have  been  described  in 
Section  IV. D.  Some  of  the  water  melted  by  the  heat  lamp  adhered  to  the 
surface  on  which  it  melted,  and  some  diffused  into  unmelted  frost  around  the 
edges  of  the  illumination.  The  test  was  not  made  quantitative  with  lamp 
intensity  calibrations.  However,  calculations  based  on  the  nominal  intensity 
values  gave  reasonably  good  estimates  for  the  time  required  for  the  frost  to 
melt.  Very  few  drops  of  water  formed  on  the  machined  SOFI.  Water  on  an  as- 
sprayed  panel  tends  to  collect  in  droplets  more  readily  than  on  a  machined 
SOFI  surface,  and  melted  frost  might  pose  more  of  a  problem  on  the  as-sprayed 
SOFI  by  refreezing  into  drops  of  ice. 

In  test  12  the  surface  temperature  of  the  SOFI  was  observed  to  decrease 
about  3.3°C/hr  (5.9°F/hr)  after  the  appearance  of  froBt.  This  3.3°C/hr  tem¬ 
perature  decrease  represents  a  1.7%/hr  change  in  the  total  temperature  drop  of 
about  190°C  across  the  SOFI,  and  therefore  in  the  heat  transfer  rate  through 
the  SOFI.  The  SOFI  surface  temperature  decreased  because  the  frost  growth 
served  as  additional  insulation  between  the  SOFI  and  the  ambient  air.  During 
the  45  min  before  frost  appeared,  Tsurf  was  approximately  25°C  lower  than 
Tamb*  Although  the  frost  surface  temperature  was  not  measured  (the  thermo¬ 
couple  was  mounted  on  the  SOFI  surface  under  the  frost),  it  would  have  been 
initially  the  same  as  the  SOFI  surface  temperature.  The  1.7X/hr  decrease  in 
the  heat  transfer  rate  through  the  SOFI  would  have  been  reflected  in  a 
1.7%/hr  decrease  in  the  25°C  temperature  difference  between  the  frost  surface 
and  the  ambient  air.  Therefore,  we  can  conclude  that  the  surface  temperature 
of  the  frost  increased  0.43°C/hr,  whereas  the  temperature  at  the  Interface  of 
the  SOFI  and  the  frost  decreased  by  3.3°C/hr.  The  temperature  drop  across  the 
frost,  therefore,  increased  by  3.7°C/hr  (6.7°F/hr)  as  the  frost  thickness  grew 
at  the  rate  of  0.063  in. /hr.  Since  the  total  heat  transfer  rate  across  the 
SOFI  was  calculated  to  be  78  Btu/ft^-hr,  the  frost  thermal  conductivity  can  be 
estimated  to  be  0.061  Btu/hr-ft-°F  or  0.11  W/m-K,  which  is  approximately  equal 
to  the  theoretical  and  experimental  values  quoted  in  Dietenberger  et  al.^  for 
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the  frost  temperature  and  density  of  test  12.  Thermal  conductivities  for  the 
other  tests  on  the  0.5-ln.  panel  have  been  calculated  and  plotted  in  Fig.  8. 
along  with  the  theoretical  values  from  Dietenberger  et  al.**  Some  uncertainty 
is  introduced  by  fluctuations  in  the  ambient  temperatures  during  the  tests. 
However,  only  one  test  point,  that  of  test  8,  lies  out  of  the  range  of  good 
agreement.  Both  the  heat  transfer  rate  and  the  frost  thickness  growth  rates 
were  smaller  for  the  tests  with  the  0.9-in.  SOFI  panel.  As  a  consequence,  the 
temperature  drop  across  the  frost  (AT  -  qd/k)  was  smaller  in  those  tests,  so 
that  the  frost  growth  did  not  perturb  greatly  the  SOFI  surface  temperature. 


31 


.f+  : 


0  0.01  0.02  0.03  0.04  0.05  0.06  0.07  0.08  0.09 

FROST  DENSITY,  g/cm3 

Thermal  conductivities  of  frost  deduced  from  tests  8  through  13 


VII.  SUMMARY 


Small-scale  experiments  have  been  performed  with  SO FI  panels  shaved  to 
0.3-  and  0.9-in.  thicknesses.  The  rear  sides  of  the  SOFI  panels  were  cooled 
to  LN2  temperatures;  the  front  sides  were  exposed  to  an  environmental  chamber 
with  controlled  temperature  and  humidity.  The  measured  rates  of  frost  forma¬ 
tion  and  the  observed  surface  temperatures  are  consistent  with  predictions 
obtained  with  a  simple  heat  balance  model.  The  measured  frost  densities,  on 

the  other  hand,  were  less  than  6  lb/ft  ,  much  lower  than  those  predicted  by 

7 

the  UDRI  frost  formation  model.  In  a  worst-case  test,  hard  ice  was  observed 
when  water  was  dripped  down  the  surface  of  a  0. 5-in. -thick,  L^-chilled  SOFI 
panel.  A  preliminary  study  of  solar  melting  produced  little  rundown  or  ice. 
Frost  densities  of  less  than  0.1  g/cm^  (6  lb/ft^)  were  obtained,  except  in  the 
test  in  which  water  was  dripped  onto  the  surface. 

Validation  by  the  present  experiments  of  predictions  based  on  a  heat 
balance  model  applies  only  to  natural  convection  in  the  laminar  flow  regime. 

A  large  fraction  of  the  flow  on  the  space  shuttle  fuel  tank  will  be  in  the 
turbulent  regime  because  of  the  height  of  the  tank.  Larger-scale  tests  would 
be  required  to  validate  the  turbulent  heat  transfer  coefficients. 
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APPENDIX.  HEAT  BALANCE  MODEL 


The  surface  temperature  Tsurj  of  the  SOFI  can  be  predicted  with  a  heat 
balance  model  In  which  <lCOIKj»  the  heat  conducted  through  the  SOFI,  la  equated 
to  the  heat  deposited  on  the  SOFI  surface  by  convective  heating  qconv*  radia¬ 
tion  qrad,  and  condensation/fusion  of  water  vapor  from  the  air  latent*  M 
follows: 


^cond  “  ^conv  +  ^latent  +  '•rad 
The  heat  conducted  through  the  SOFI  is  given  by 

^cond  "  *2  "AX  ^surf  ”  Tcryo^  Btu/ft  -hr 

where 


(A-l) 


(A-2) 


AX  »  SOFI  thickness  in. 


avg 


0.0017  +  2.45  x  10' 


^Tsurf  +  Tcryo^ 


,-5 


avg 


Btu/hr-f t-°R 


For  example,  with  Tsurf  -  49l°R  (26°F),  Tcryo  -  138°R,  and 
AX  -  1  in. ,  then 


k  -  0.0094  Btu/hr-f t-°R 


and 


9Cond  "  40  Btu/ft2-hr 
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The  equation  for  convective  heating  la  given  by 


<lconv  “  hconv  AT  Btu/ft2-hr 


(A-3) 


where 


“  ^amb  “  ^surf 


For  laminar  convection  over  a  vertical  wall,^ 


hconv  “  °-29  (^)0,25  Btu/f t2-hr 


L  -  length  of  wall  ft 


hconv  *  °*22  (AT)0,25  Btu/f t2-hr  for  L  -  3  ft 


The  equation  for  radiation  la  denoted  aa 


«rad  "  *E  “  Tsurf>  Btu/ft2-hr 


(A-4) 


where 


E  -  emissivity  *  1 


o  *  Stephan-Boltzmann  constant  ■  1.714  x  10~®  Btu/ft2-hr-°R^ 


NOTE:  q  can  be  expanded  in  terms  of  AT  to  give 


"rad  "  E  TLh  1T  (1  '  31T/2T„b  •••> 
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The  equation  for  water-vapor  condensation/fusion  is  given  by 


where 


^latent 


h 

cony 

C  4 
P*air 


<Wamb  “  Wsurf>  hfg 


Btu/f t2-hr 


(A-5) 


hconv  “  convective  heat  transfer  coefficient 

Cp  a^r  “  specific  heat  of  air  at  constant  pressure 
-  0.24  Btu/lb-°F 

h^  ■  heat  of  sublimation 

■  675  cal/g  m  1220  Btu/lb 

Wamb  *  am^ient  water-vapor  concentration,  lb  l^O/lb  air 

wsurf  “  saturated  water-vapor  concentration 

at  surface  temperature,  lb  t^O/lb  air 

W  can  be  related  to  (Pamb  -  Psurf)  (the  partial  pressure  of  water,  in  psl)  by 

W  -  0.042P 

Therefore 

Olatent  '  *7  (AT)0'25  (Pamb  -  P8urf>  Btu/ft2-hr 
The  mass  rate  of  frost  accumulation  is  written 

m  -  0.0388  (AT)0*25  (PaBb  -  Psurf)  lb/ft2-hr 

-  lb/ft2-hr  (A-6) 
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Substituting  the  above  expressions  into  Eq.  (A-l),  one  obtains  the  heat 
balance  equation,  with  AT  as  the  unknown.  Although  the  equation  is  nonlinear, 
a  rapid  convergence  is  obtained  with  an  iterative  procedure.  A  computer  pro¬ 
gram  written  in  BASIC  and  a  sample  calculation  are  included  here  (Pig.  A-l). 
The  computer  program  is  written  with  the  temperatures  in  degrees  Celsius  (°C) 
instead  of  degrees  Rankine  (°R). 


\ 

F 

i 


i 


»LIST 


10  REN  PROGRAM  CALLED  T SURF-2  REVISED  15/ JUNE/82 

20  INPUT"PAMB='\  PI;  :  PRINT'‘PSI 
2©  INPUT"  TAM6=">  Ti;  :  PRINT"DEG.  C"; 

25  INPUT"  EMISSIVITV=%  El; 

40  INPUT"  S0FI  THICKNESS=“,D2 

50  C1=0.  46 

60  C5=29/D2 

70  B=Tl+272 

80  C2=Ei*7.  2E-08*6**2 

90  C4=2/<2*B:> 

100  D2=20  :  REM  FIRST  GUESS  FOR  TAMB=TSURF 
120  FOR  1=1  TO  20 

120  T4=Tl-D2+272  :  REM  SURFFiCE  TEMP  CfiLC 

140  P2=0.  08S6*EXP<22.  46-6121/T4> 

150  C6=D2**0.  25 

155  C2=54.  4*<P1-P2> 

160  D1®KC5-C2*C6>/<C1*C6+C2*<1-C4*D2>  > 

170  IF  ABS<D1-D2><0.  01  THEN  180 

172  IF  1=20  THEN  270 

174  02=D1 

176  NEXT  I 

18©  PRINT  .  PRINT 

181  PRINT"  PfiMB  =  ";P1;"  PSI  PSURF  =  P2j 

182  PRINT"  PSI  TAMB  =  "i  Tl; "  DEGREES  C" 

182  T2=T1-D1 

184  D=0.  0219*D1**0  25*<P1-P2> 

190  PRINT"  TSURF="i T2; "  DEGREES  C  DM/DT  =  "i  D 
200  0 

210  PRINT"  TfiMB-TSURF=  ";  D1 

220  <?:<?:  0 
220  C« 

260  STOP 

270  PRINT  "FAILED  TO  CONVERGE  LfiST  VALUE  =  ";D2;,,NEW  VALUE*"; 
280  END 


Fig.  A-l.  BASIC  computer  program  for  surface  temperature. 
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TEST  12 


»PR  I  NT  "SAMPLE  CASE 
SAMPLE  CASE  TEST  12 

»RUN 

PAM6=0.  2217PSI  TAM6=14D£G.  C  EMISSIVITV=1.  6  SOFI  THICKNESS*©.  5 


PAM6  =  B.  2217  PS1  PSURF  =  0.  042911878888199  PS  I  TAMB  *  14  DEGREES  C 
TSURF=-8.  290808792211  DEGREES  C  DN/DT  =  8.  9280202928777E-02 

TAMB-TSURF*  22.  290808792211 


*>*<*260  STOP*** 


Fig.  A-l.  BASIC  computer  program  for  surface  temperature.  (Cont.) 
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The  Ukonuqr  OftntloM  of  Thu  teraptu  Corporation  ta  conducting  aopor- 
lantil  and  theoretical  inutl|itlHi  hchhi;  for  the  ralwilM  tad  applica¬ 
tion  of  oclaatiflc  advance*  to  aao  nllltary  apaca  ayataaa.  Voraatlllty  aad 
flexibility  han  kaoa  developed  to  a  high  dagrac  bp  tba  laboratory  poroooaol  la 
daallag  ultb  tba  aaay  problaa*  aacoaatarad  la  tba  aatloa'a  rapidly  developing 
apaca  ayataaa.  Expart lac  la  tba  lataat  aclaatlfic  dovalopaant*  la  vital  to  tba 
accoopllehnent  of  taaba  ralatad  to  tboaa  problene.  Tba  laboratorlaa  that  con- 
trlbuta  to  tbla  raaaareb  ara: 

Aaroahralca  laboratory;  Lauach  ooblcla  aad  raoatry  aarodynaalca  and  haat 
tranafar,  ptopulalon  cKSXStry  aad  fluid  aachaalca,  atructural  aaebanlca,  flight 
dyaaalca;  blgb-taoparatura  t  ha  ran aachaalca ,  gaa  klaatlea  aad  radlatloa;  roaoarch 
la  anvlroaaantal  cbaalatry  aad  eeotMlnatloa;  cu  aad  pul  aad  chaalcal  la  oar 
daaalopaaat  Including  chaalcal  klaatlea,  apactroaeopy,  optical  raaoaatora  and 
baa  pointing,  ataoapbarlc  propagation,  laaar  affacta  aad  countoiaaaaure*. 

Chaulatry  aad  Phyalca  Laboratory;  ataoapbarlc  chaalcal  roactloaa,  atao- 
apbarlc  optica,  light  acattarlng,  atata-apaclflc  chaalcal  roactloaa  aad  radia¬ 
tion  traaaport  ta  tackat  ploaaa,  applied  laaar  apactroaeopy,  laaar  cbaalatry, 
battery  elactrockanlatry,  apaca  vacuua  aad  radlatloa  affocta  on  aaterlala,  lu¬ 
brication  aad  aurfaca  phauonena,  tberaloalc  aalaaloa,  pbotoeonoltlvc  aaterlala 
aad  datactera.  ataalc  froguaacy  ataadarda,  aad  bloaavlraaaaatal  raaaareb  and 
aaaltorlag* 

gfoctrealco  baoaarcb  Laboratory  i  mcroalactroalca.  Coda  lov-nola*  aad 
poaar  Juvlcoo,  aoalcooductor  laeaTa,  alactraaagaatlc  aad  optical  propagation 
phaaoaaaa,  paata  alactroalca,  laaar  coanaaleatiaae,  Ildar,  aad  olactro-optlca; 
coaaaalcat lea  acloacaa,  applied  alactroalca,  aaalcooductor  cryetal  aad  device 
p by a lea,  radlooatrlc  laaglag;  allllaotor  oavo  aad  alcroaaoa  technology • 

luforaatloa  bclaacaa  baaaarcb  Offltoi  Pragraa  varlflcatloa,  pragraa  trana- 
latlauj  porforaaaco-aaaa It loo  ayataa  Jdaiga,  dlatrlbutad  archl toctorao  for 
apocaharaa  oaoputara,  fault-tolar aat  cooputar  ayataaa,  artificial  latalllgaaca, 
aad  alcraalactroalca  applleatlaaa. 

Material*  bclaacaa  laboratory:  Braaloyaaat  of  non  aaterlala:  natal  aatrlu 
cnapmltoo.  polyaara.  aaa  aaa  feme  of  carbaa;  caapoaoat  failure  aaalyala  aad 
rallabllltyi  fractara  aachaalca  aad  atraoe  corroaloa;  aualuatlon  of  aatarlala  la 
apaca  eavlreMcati  aaterlala  parfomaaca  In  apaca  troaapartatloa  ayataaa;  aaal- 
yola  of  ayataaa  aalaarabllity  aad  aurvlvablllty  la  allay  Induced  aaolraaaanta. 

beat*  gcloacoa  laboratory:  Ataoapbarlc  aad  laaaapbarlc  phyalca,  radlatloa 
f ran  tba  ataoapbara,  aanaity  aad  coapoeltlea  af  tba  upper  otaoopbaro,  aurora* 
aad  alrglaat  nagaataaphorte  phyalca,  caaalc  raya,  goaaratlaa  aad  propagation  of 
ploaaa  *****  la  tba  aagaotoopbaroj  aolar  phyalca,  lafrarod  aatraaaap;  tba 
affocta  af  naalaar  aaploalaaa,  aagnattc  atoraa,  aad  aolar  actlalty  an  tba 
earth’*  ataoapbara,  laaaaphara,  and  aagaataaphara;  tba  affocta  at  aptlcal, 
alactraaagaatlc,  aad  particulate  radlatlao*  la  apaca  aa  apaca  ayataaa. 


